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We compute the cross section and differential distributions for the production of a Z boson in asso¬ 
ciation with a hadronic jet to next-to-next-to-leading order (NNLO) in perturbative QCD, including 
the leptonic decay of the Z boson. We present numerical results for the transverse momentum and 
rapidity distributions of both the Z boson and the associated jet at the LHC. We find that the 
NNLO corrections increase the next-to-leading order (NLO) predictions by approximately 1% and 
significantly reduce the scale variation uncertainty. 


PACS numbers: 12.38Bx 

The Drell-Yan production of lepton pairs is a bench¬ 
mark process at hadron colliders like the Large Hadron 
Collider (LHC). The production of Z bosons (or off-shell 
photons) with subsequent leptonic decays has both a 
clean and readily identifiable signature and a large event 
rate. It is a key process for precision measurements of 
electroweak (EW) parameters, and also allows to probe 
various aspects of the strong interaction, including par- 
ton distribution functions (PDFs), the strong coupling 
constant and the behaviour of processes involving 
multiple scales. It is moreover a key ingredient in cal¬ 
ibrating several parts of the detector (including the jet 
energy scale) and can potentially be used to measure the 
luminosity of the collider. At the LHC, the Z boson is 
almost always produced together with additional QCD 
radiation thereby providing a perfect testing ground for 
our theoretical understanding of both strong and elec¬ 
troweak physics in a hadronic environment. Together 
the combination of precise experimental data and reli¬ 
able theoretical predictions enables a variety of precision 
measurements at the LHC mm- 

The importance of the neutral current Drell-Yan pro¬ 
cess is also reflected in the effort to make the theoretical 
predictions as precise as possible. For inclusive Z pro¬ 
duction, theoretical predictions at per-cent level accuracy 
are available. To attain this level of precision, a variety 
of higher-order corrections in QCD and the EW theory 
had to be considered. The cross section for Z produc¬ 
tion is known at next-to-next-to-leading order accuracy 
(i.e. at two loops) with respect to QCD corrections [3]. 
Corrections beyond this order have been studied in the 
soft-virtual approximation [4]. The NNLO QCD correc¬ 
tions have been combined with a resummation of next-to- 
next-to-leading logarithmic effects [5] which is necessary 
to predict the transverse momentum distribution of the Z 
boson at small Pt and matched with parton showers [ 6 ] . 
In the electroweak theory, the next-to-leading order cor¬ 
rections m E] and the mixed QCD-EW corrections [9] 
also contribute to the precise description of this process. 



FIG. 1. Representative Feynman diagrams for (a) two-loop 
Z boson-plus-three-parton amplitudes (b) one-loop Z boson- 
plus-four-parton amplitudes and (c) tree-level Z boson-plus- 
five-parton amplitudes. 

Drell-Yan production in association with hadronic jets 
has also been intensively studied. The NLO QCD cor¬ 
rections for Z + 1 jet m, Z + 2 jets m, Z + 3 jets m 
and Z -b 4 jets m are known while the NLO EW cor¬ 
rections for Z + 1 jet m and Z + 2 jets m have also 
been derived. 

In this letter, we report on the calculation of the NNLO 
contributions to the neutral-current Drell-Yan process in 
which the dilepton pair is produced in association with a 
hard, visible hadronic jet, 

pp Zj^'^ + jet ^ + jet + X. 

Our results are obtained in the form of a parton-level 
event generator that provides the corrections in a fully 
differential form, including the Z/ 7 * boson decay to two 
charged leptons. The final state of the hard-scattering 
process is completely reconstructable and the application 
of an invariant mass cut on the lepton pair can ensure 
that the process is dominated by resonant Z bosons. 

The NNLO corrections to Z boson + jet production 
in hadronic collisions receive contributions from three 
types of parton-level processes: (a) the two-loop cor¬ 
rections to Z boson-plus-three-parton processes [16], (b) 
the one-loop corrections to Z boson-plus-four-parton pro¬ 
cesses El [m and (c) the tree-level Z boson-plus-five- 
parton processes pn [TS] . Figure shows representative 
Feynman diagrams for each of the partonic multiplici¬ 
ties. The ultraviolet renormalized matrix elements for 
these processes are integrated over the final state phase 
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space appropriate to Z boson + jet final states. All three 
types of contributions are infrared-divergent and only 
their sum is finite. While infrared divergences from the 
virtual corrections are explicit in the one-loop and two- 
loop matrix elements, divergences from unresolved real 
radiation become explicit only after phase space integra¬ 
tion. The divergences are regulated using dimensional 
regularization, and a variety of methods have been used 
for their extraction from the real radiation contributions. 
All these methods are based on the isolation of the diver¬ 
gent configurations, which are then integrated over the 
phase space and added to the virtual corrections to yield 
a finite result: sector decomposition m, sector-improved 
residue subtraction [20], antenna subtraction [21], qt- 
subtraction [22] and N-jettiness subtraction [23] have all 
been applied successfully in the calculation of NNLO cor¬ 
rections for a range of LHC processes. 

In this calculation we employ the antenna subtraction 
method [21] in which the real radiation subtraction terms 
are constructed from antenna functions. These antenna 
functions capture all the unresolved radiation emitted be¬ 
tween a pair of hard radiator partons. For hadron-collider 
observables, either hard radiator can be in the initial or 
final state, and all unintegrated and integrated antenna 
functions were derived in Refs. [25H28] . The cross section 
corresponding to an initial partonic state ij is given by, 

[^cf^%nlo ~ ^^ij,NNLo\ 

\_^^^^NNLO ~ ^^Ij,NNLo\ 

[^^Yj^NNLO ~ ^^Yj,NNLo\ 5 ( 1 ) 

where each of the square brackets is finite and well be¬ 
haved in the infrared singular regions. The construc¬ 
tion of the subtraction terms follows closely 

the procedure established for jet production [29] and 
Higgs + jet production [30]. Powerful checks of our for¬ 
malism are that (a) the poles in the dimensional regular¬ 
ization parameter e cancel analytically and (b) that the 
subtraction terms accurately reproduce the singularity 
structure of the real radiation matrix elements. 

Using the antenna subtraction method, we have de¬ 
rived the corresponding subtraction terms for all partonic 
initial states and all color contributions for Z boson-plus- 
jet production through to NNLO and implemented them 
in a parton-level event generator. With this program, we 
can compute any infrared safe observable related to Z + 
jet final states to NNLO accuracy. The Z boson decay to 
two charged leptons is included, such that realistic event 
selection cuts on the leptonic final state can be applied. 
Renormalization and factorization scales can be chosen 
(dynamically) on an event-by-event basis. 

For our numerical computations, we use the 
NNPDF2.3 parton distribution functions m with the 
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corresponding value of as{Mz) = 0.118 at NNLO, and 
Mz = 91.1876 GeV. Note that we systematically use 
the same set of PDFs and the same value of as{Mz) 
for the LO, NLO and NNLO predictions. The fac¬ 
torization and renormalization scales are chosen to be 
ji = jiF = jiji = Mz, with a theoretical uncertainty esti¬ 
mated by varying the scale choice by a factor in the range 
[ 1 / 2 , 2 ], 

We require that the leptons have pseudorapidity, \r]^\ < 
5 and that the dilepton invariant mass is close to the Z 
boson mass, 80 GeV < m££ < 100 GeV. Jets are recon¬ 
structed using the anti-/cT algorithm [32] with R = 0.5 
and are required to have > 30 GeV and < 3. 
With these cuts, we find that the total cross section at 
different perturbative orders is given by, 

o-LO = 103.6^11 pb , 
c^NLO = 144. pb , 

O'NNLO = 145.8l5‘.2 (2) 


so that the inclusive NNLO corrections amount to a 1% 
increase on the NLO cross section. 

More information on the impact of the NNLO QCD 
corrections can be gained from differential distributions 
in the kinematical variables of the Z boson and the jet. 
In the kinematical distributions and ratio plots, the error 
band describes the scale variation envelope as described 
above, where the denominator in the ratio plots is eval¬ 
uated at fixed central scale, such that the band only re¬ 
flects the variation of the numerator. Figure [^a) shows 
the inclusive leading jet transverse energy distribution in 
10 GeV bins. Due to the inclusiveness of the observable, 
events with two or three jets with > 30 GeV and 
< 3 are also included. The relative corrections are 
further exposed in Figure j^b) where we show the ratio, 
K = {p = M^). The band shows 

the effect of varying p in the range [1/2, 2]Mz in the nu¬ 
merator while keeping p = Mz in the denominator. For 
our set of cuts and input parameters, we see that the 
NLO corrections increase the cross section by between 
30% to 70%. At low transverse momentum the NNLO 
corrections are a positive correction of approximately 1%. 
The variation with the unphysical scales is significantly 
reduced as we move from NLO to NNLO. 

The rapidity distribution of the leading jet is displayed 
in Figure Note that the distribution is restricted by 
the requirement that < 3. We see that the NLO cor¬ 
rections are typically 35%-40% and relatively fiat. The 
NNLO corrections increase the cross section by approxi¬ 
mately 1% over the whole range of y^^^ with a significantly 
reduced scale dependence. 

The Z boson pt distribution in inclusive Z + jet pro¬ 
duction is shown in Figure!^ where we observe an inter¬ 
esting structure around p^^ 30 GeV. This behaviour 
arises from the fact that the Z boson is recoiling against 
a complicated hadronic final state that contains at least 
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FIG. 2. (a) Transverse momentum distribution of the lead¬ 
ing jet in inclusive Z + 1 jet production in pp collisions with 
= 8 TeV at LO (blue), NLO (green), NNLO (red) and (b) 
Ratios of different perturbative orders, NLO/LO (turquoise) 
and NNLO/NLO (mauve). 



yjet 


FIG. 3. (a) Rapidity distribution of the leading jet in inclu¬ 

sive Z + 1 jet production in pp collisions with ^/s = 8 TeV 
at LO (blue), NLO (green), NNLO (red) and (b) Ratios 
of different perturbative orders, NLO/LO (turquoise) and 
NNLO/NLO (mauve). 


one jet with > 30 GeV. For this set of cuts, the lead¬ 
ing order process is constrained to have pf > 30 GeV, 
while higher order real radiation corrections lift this limi¬ 
tation, since extra partonic radiation can also balance the 
transverse momentum of the leading jet. This Sudakov 
shoulder phenomenon is also observed in H+jet produc¬ 
tion [Milan]; it is well understood [33] and leads to large 
higher order corrections, which require logarithmic re¬ 
summation. Nevertheless, the NNLO corrections tend to 
stabilise the NLO result, and in fact simply represent a 
NLO correction to the p^ distribution for Z + jet events 
in this region. At larger transverse momenta, the NNLO 
corrections increase the prediction by approximately 1%. 

Figurej^shows the rapidity distribution of the Z boson. 
The NLO and NNLO corrections are largest in the for¬ 
ward/backward regions where the phase space is enlarged 
by the possibility that the hadronic radiation partially 
balances leading to a smaller Z In these regions, one 
of the parton momentum fractions is reaching a maximal 
value. In the central region, the NNLO corrections are 
very small with a reduced scale dependence. 

In the differential distributions we observe that the cor¬ 
rections are not always uniform, implying that a rescaling 
of lower-order predictions is insufficient for precision ap¬ 
plications. The need for using the fully differential higher 
order predictions can be understood for example in the 
extraction of parton distributions functions from Z + jet 
production. At leading order, the momentum fractions 


of the incoming partons is completely fixed by the trans¬ 
verse momenta and rapidities of the final state particles. 
At higher orders, the real radiation spoils the leading 
order kinematics, such that 

V( 

(^y/p/V+TO^ exp{-yz) + p^T exp(-2/jet)^ , 

where the equality is restored only for the leading order 
kinematics {p^ = The relevant x ranges probed 

by Z boson-plus-jet production is thus determined by the 
transverse momentum and rapidity distribution of the Z 
boson and the jet. For our cuts, the smallest momentum 
fractions probed are x ^ 8 -10 “^, and smaller values of x 
can be attained by enlarging the rapidity interval or by 
lowering the transverse momentum cut. 

In this manuscript we have presented the complete 
NNLO QCD calculation of Z boson production in as¬ 
sociation with a jet in hadronic collisions including all 
partonic subprocesses. This process is measured exper¬ 
imentally to high precision mm and is an important 
ingredient to a variety of precision studies of Standard 
Model parameters and derived quantities as well as a 
key element in the LHC detector calibration. We have 
achieved this using the antenna subtraction method that 
has been successfully applied to other processes at the 
LHC. For all of the observables considered here, we ob- 


exp( yz)+p’T ®xp( 
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FIG. 4. (a) Transverse momentum distribution of the Z bo¬ 
son in inclusive Z + 1 jet production in pp collisions with 
= 8 TeV at LO (blue), NLO (green), NNLO (red) and (b) 
Ratios of different perturbative orders, NLO/LO (turquoise) 
and NNLO/NLO (mauve). 


served a very significant reduction of the respective un¬ 
certainties in the theory prediction due to variations of 
the factorization and renormalization scales with a resid¬ 
ual NNLO scale uncertainty of around 1% on the nor¬ 
malization of the distributions. Our calculation will be a 
crucial tool for precision studies of Z boson -b jet produc¬ 
tion in the upcoming data taking periods at the CERN 
LHC. 

This research was supported in part by the Swiss Na¬ 
tional Science Foundation (SNF) under contracts 200020- 
149517 and CRSII2-141847, in part by the UK Sci¬ 
ence and Technology Facilities Council as well as by 
the Research Executive Agency (REA) of the Euro¬ 
pean Union under the Grant Agreement PITN-GA-2012- 
316704 (“HiggsTools”), and the ERC Advanced Grant 
MC@NNLO (340983). 

Note added: After this paper was initially submit¬ 
ted to Physical Review Letters, a second calculation (em¬ 
ploying a different subtraction scheme) of Z+jet produc¬ 
tion at NNLO precision has been presented and pub¬ 
lished [34]. In coordination with the authors of [34|, 
we performed an in-depth comparison, by running our 
code with their settings (cuts, parton distributions, scale 
choice). This comparison uncovered an error in the nu¬ 
merical code used in [34] , which alters their published re¬ 
sults. After correction of this error, the code developed 
in [34] agrees with our results. 



FIG. 5. (a) Rapidity distribution of the Z boson in inclusive 
Z + 1 jet production in pp collisions with ^/s = 8 TeV at LO 
(blue), NLO (green), NNLO (red) and (b) Ratios of different 
perturbative orders, NLO/LO (turquoise) and NNLO/NLO 
(mauve). 
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